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1. INTRODUCTION 


The variations of the upper atmosphere are due to a wide variety of 
sources and cover a broad range of frequencies and size scales. The shortest 
period, smallest scale variations are due to turbulence. Certain variations 
of only slightly longer period and larger scale are thought to be due to the 
action of gravity waves. For a review of the evidence for gravity waves and 
the possible mechanisms for generation and dissipation of gravity waves see 
the previous report by Justus and Woodrum^^ . There has been no unambiguous 
resolution of gravity waves in the upper atmosphere (i.e. simultaneous obser- 
vation of amplitudes, phases, and frequency sufficient to verify propagation 
according to the theoretical dispersion equation). However, there is strong 
circumstantial evidence for their occurrence and the short period (i.e. <24 
hr) irregular variations discussed in this report may be assumed to result 
from gravity waves. 

Atmospheric tides cause variations with periods of 24 hours and har- 
monics (e.g. 12 and 8 hours). The strongest tidal influence comes from 
solar heating, but lunar tides also exist similar to those produced in the 

ocean. For a review of theory and observations of atmospheric tides see the 

( 2 ) 

book by Lindzen and Chapman 

Longer period phenomena are the result of synoptic variations and plane- 
tary scale waves. Although planetary waves have been resolved up to about 
(3) 

30 km , their propagation to higher altitudes has been somewhat uncertain 
prior to this report. Still longer period phenomena are the seasonal, semi- 
annual, annual, and quasi-biennial oscillations. 



This report is a continuation and revision of results reported earlier^. 

Two analysis techniques are used to study gravity wave, planetary wave, and 

long period variations* The first method, known as* the daily difference tech- 

(4 5) 

nique was developed * for application to upper atmosphere data where 
specific removal of long period and tidal variations would be difficult be- 
cause of lack of data. This approach allows for the estimation of magni- 
tudes, probability distributions, and vertical structure functions of the 
irregular variations of less than one pay period. It also can be used to 
determine the magnitude of the variation* due to planetary waves of periods 
up to a few days. The second data analysis technique is harmonic analysis 
of time series data. This approach can be used, when sufficient data exists, 
to estimate tides, and long period variations such as annual and quasi- 
biennial oscillations. 

The purpose of this report is two fold: first to present the results of 
an investigation of the effects of gravity waves, tides, and planetary waves 
in producing upper atmospheric variations, and second to summarize the re- 
sults in such a fashion that they can be used in engineering design problems, 
in particular the effects of the atmospheric variations on the stability, 
control, and dynamic heating of the space shuttle vehicle. 
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2. DATA ANALYSIS AND INTERPRETATION 


Structure Functions . An alternate form of the correlation function known 
as the structure function was first used extensively by Russian meteorologists 
in the analysis of turbulence. The structure function of a statistically stat- 
ionary time varying process f(t) is given by 

D(t) = <[f(t + x) - f (t) ] 2 > (1) 

where the structure function D depends only on the time displacement x because 
of the statistical stationarity . The angle brackets in (1) denote averaging. 

Structure function analysis can be applied to wave phenomena also. Con- 
sider the process f(t) to be a wave of frequency co and amplitude A, i.e. f(t) = 

A sin (cot + a) , where a is a phase angle. If the averaging process is taken 

to be integration over any whole number of periods, where the period T * 27r/co, 

2 2 

then the mean square value of f(t) is <f > = A /2. It is also easy to deter- 

2 

mine that the structure function in this case is given by D(t) = A (1 - cos cox) 

Thus the structure function of a single component wave field has the properties 

2 

of being proportional to x for very small time displacement x since the lead- 
ing term in 1 - cos cox is co^x^/2, and being periodic with the same frequency co 
as the wave. Notice that the first maximum in D(x) would occur at i = tt/co * T/2 
so the first maximum occurs at the half period value. The value of the maximum 
is 2A 2 . 

The structure function analysis of waves can easily be extended to a set 
of waves which are part of a harmonic series. Consider f(t) to be periodic 
with period T and to be made up of N different waves each wave n having a fre- 
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quency = 2mr/T, i.e. 


f(t) «. Z A sin (w t + a ) (2) 

n n n 

where the summation is from n = 1 to n = N. Using the same type of averaging 
by integrating as on the single wave, one can determine that the mean square 
value of f is given by 


2 2 
<f > = Z A /2 

n 


(3) 


and that the structure function is given by 


D(x) 



(1 - COS 0) t) 

n 


(4) 


Thus again the behavior at small time displacement is proportional to x , since 

2 2 2 

the leading term in (4) is Z A a) x /2. The structure function in this case 

n n 

is again periodic with period T and has components with amplitudes related to 
the square of the amplitudes of the corresponding frequencies in the process 
f(t). 

If the process f(t) is made up of sufficient components that it must be 
considered as having a continuous spectrum of amplitudes A(to) (as is the case 
with turbulence and may also be the case when gravity wave modes have signi- 
ficant nonlinear modal interaction) then f(t) can be represented by 


f(t) 


. , N iwt , 
A(oj) e do> 


(5) 


(For a more rigorous discussion of spectral representation see Lumley and 
Panofsky^^ and Lin^\ The mean square value of f(t) is given by 


4 



<f^> = 2 | 4>(w) da) (6) 

where <(> (cu) is the spectral density of the mean square variations of f(t). The 
structure function is given by 

D(x) = 2 

o 

Again the initial variation of D(t) for small time displacements is proportional 

_ 2 . 
to t since 

D(t) - [ 

o 

However, at large time displacements D(x) is not periodic but approaches a con- 

2 

stant value of 2<f >. This is seen readily from the fact that D(t) is related 
to the correlation function p(x) by 

D(t) = 2<f 2 > [1 - p (t) ] (9) 

for all stationary processes f(t). Hence, because p(x) approaches 0 for large 
x , then 

D(t) = 2 <f 2 > t ^ » (10) 

In addition, if there is an appreciable range of frequencies over which (co) 
varies as a power law 

(J> ((jo) = C 0) a (11) 

then there is a range of time displacements over which 


r 


2 2 
u <Kw) doi] x 


( 8 ) 




<f>(w) (1 - cos ojt) dto 


(7) 
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D(t) ~ C[ x a (1 cos x) dx] x a 1 (12) 

o 

where x = cox and the integral in brackets is some definite constant value. 
Thus observed power law behavior of the structure function (i.e. a - 1) can 
be related to the power law of the spectral density (i.e. -a). 

As an example of how the integral time scale T = 

o 

tained from the time structure function, consider the case in which the 
correlation p(x) of the process f(t) is given by exp (- x/T). For this par- 
ticular case the small time behavior of the structure function is D(x) ~ 

2 

2<f > x/T. If this early time period slope is extrapolated to intersect 

2 

with a similar extrapolation of the long time variation, D(x) ~ 2<f >, the 
graphical intersection point occurs at a value x = T. Similar graphical 
techniques can be devised to determine integral scales if different corre- 
lation functions, with different early time behavior, apply. Also similar 
techniques can be used with spatial structure functions to estimate size 
scales (see section 8). 

Daily Difference Analysis . The daily difference analysis technique was 
developed^’*^ for applications where limited data did not allow explicit 
separation of the tidal components in order to determine the small scale 
irregular variations. As an example of the application of this technique 
consider a vertical profile of a parameter F(z, t) over height z at time t 
where F may be a wind component, pressure, density, or temperature. We con- 
sider that F is made up of a prevailing value F q (z ) which is time invariant, 
plus a long period (e.g. seasonal, annual, or quasi-biennial oscillation) 
component S(z, t), a planetary scale or synoptically varying component 
P(z, t), a tidal component T(z, t), a gravity wave or short period irregular 


p(x) dx can be ob- 
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component G(z, t) , and a still smaller scale component made up of measurement 
error and turbulence E(z, t) . Thus 

F(z, t) = F q (z) + S(z, t) + P(z, t) + T(z, t) 

+ G(z, t) + E(z, t) (13) 

The component P(z, t) would be composed of traveling waves only, all truly 
standing waves would be included in the component F q (z) or seasonally fluct- 
uating standing waves would be included in S(z, t). We now choose two pro- 
files of F at times t^ and t 2 such that t 2 - t^ = At = 24n hours where n is 
an integer. If, at any selected altitude z, we difference the corresponding 
value of F then 

AF n (z) = F(z, t 2 ) - F(z, t 1 ) = [S (z, t 2 ) - S (z , t^) ) 

+ [P(z, t 2 ) - P(z, t x )] + [T(z, t 2 ) - T(z, t x )] 

+ [G(z , t 2 ) - G(z, t x )] + [E(z, t 2 ) - E(z , tj)] (14) 

We now make the following assumptions: (1) Assume that n is sufficiently small 
that S(z, t ? ) - S(z, t^) (i.e* n is a small number of days compared to times 
over which appreciable seasonal variation would occur) . In the analysis we 
restricted n to 15 or less days* (2) Assume that because the tidal component 
is diumally repeating and At is a multiple of 24 hours that T(z, t * 

T(z, t^) . (Any systeramatic or synoptic variation in the tidal parameters 
would be included in the component P and the seasonal variation of the tides 
would be included in the component S) . (3) The planetary scale, gravity 
wave and error components are uncorrelated with each other and are corre- 
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lated only with themselves (autocorrelation). Equation (13) can now be 
squared and averaged over an ensemble of different profile pairs all hav- 
ing the same time separation At. The result is 

<[AF n (z)] 2 > - <[P(z, t 2 ) - P(z, t 1 )] 2 > 

+ <[G(z, t 2 ) - G(z, t x )] 2> + <[E(z, t 2 ) - E(z, t 1 )] 2 > (15) 

The cross product terms in (15) have dropped out because of assumption 3 
above. If equation (15) is now expanded and the mean square values of P, 

G, and E are assumed to be independent of time (i.e. statistically station- 
ary) then the mean square data differences become 

<[AF n (z)] 2 > = 2<P 2 (z) > [1 - p p (At) ] 

+ 2<G 2 (z)> [1 - pg(At) ] + 2<E 2 (z)> [1 - Pp(At) ] (16) 

where p , p , and p are the time autocorrelation functions of P, G, and E 

r Vjr Ei 

respectively. The following assumptions are now made about the periods of 
the various remaining components: (1) the gravity wave, error and turbu- 
lence components are uncorrelated for all time differences of 1 day or 
more (i.e. n > 1), (2) the planetary wave component is of such a long per- 
iod that pp (At) - 1 for At * 1 day, but for large n the planetary wave 
component also becomes uncorrelated. Thus for single day differences, 
equation (16) becomes 

<[AF 1 (z)] 2 > = 2[<G 2 (z)> + <E 2 (z)>] (17) 
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that is, the mean square differences in the observed data are equal to twice 
the mean square magnitude of the gravity wave component (plus any contribut- 
ion from measurement errors or small scale turbulence) . For time separat- 
ions of many days (n large, say approaching 15) and under the above assum- 
ptions equation (16) becomes 

<[AF n (z)] 2 > = 2[<P 2 (z)> + <G 2 (z)> + <E 2 (z)>] (18) 

thus, at longer time separations the magnitude of the planetary wave contri- 
butions becomes added. At intermediate time separations progressively 
larger portions of the planetary wave contribution (through the factor 
1 - Pp(At)) become added. Equations (17) and (18) can be subtracted, which 
yields 


<[AF (z)] 2 > - <[AF 1 (z) ] 2 > = 2<P 2 (z) > (19) 

n l 

This allows an estimate of the contribution of planetary waves directly from 

the observed daily differences of measured data and the estimate is unbiased 

2 

with respect to the error component <E > since that component cancels in the 
subtraction process. Note, however, that this method, like any single site 
method, does not resolve the standing planetary wave components, only the 
traveling components. 

The assumptions outlined above regarding relative periods of the grav- 
ity wave and planetary wave components (and the implicit assumption that 
the errors are sufficiently small that meaningful results can be obtained from 
the analysis) are subject to verification. The results presented later in 


9 



I 


this report do, indeed, confirm these assumptions. 

The vertical structure function of the daily differences can also be 
related to the vertical structure function of the gravity wave component. 
The vertical structure function of the 1 day differences in the data 
values is given by 


D AF 1 ( ° = <UFj.Cz> C) - AF 1 (z)] 2 > (20) 

Subtraction of the expression 

AF^z) = [G(z, t 2 ) - G(z, tj.)] + [E(z , t 2 > - E(z, tj.) ] (21) 

and the comparable expression for AF^(z + O yields, after some rearrange- 
ment 

D af (?) = <[G(z + ?, t 2 ) - G(z , t 2 )] 2 > 

+ <[G(z + £ , t^) - G(z, t^)] 2> + <[E(z + £, t 2 ) - E(z, t 2 )] 2 > 

+ < [E(z + 5 , tj.) - E(z, t x )] 2 > (22) 

Several cross product terms in (22) have dropped out because of the lack of 
correlation between G and E and the lack of time autocorrelation of G and E 
over periods of 1 day. If the G and E fields are vertically homogeneous, 
then (22) will depend on £ only, otherwise it will depend separately on z 
and 5 . Under the previously employed assumption of statistical stationarity , 
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and under the assumption that the E values are uncorrelated over the height 
separation ?, (22) becomes 

(?) = 2<[G(z + ?) - G(z) ] 2 > + 4<E 2 > 

AF 1 

= 2 D g (?) + 4<E 2 > (23) 

Thus the vertical structure function of the single day differences in the 
data profiles is twice the vertical structure function of the gravity wave 
component (plus a contribution from error and small scale turbulence) . 

A similar analysis for vertical structure functions of daily differences 
of data profiles for At a large number of days, (but small compared to the 
periods in the component S) shows that 

d af (c) = 2D P ( ° + 2 d g ( c) + 4<e2> (24) 

n 

The structure function of the planetary wave component can be obtained by 
differencing (24) and (23)* 


Q^Ui-D ( « - 2D P U) (25) 

n 1 

In addition to vertical structure functions, horizontal structure functions 
of the gravity wave components can also be calculated by differencing data from 
the up and down trail trajectories of rockets or from closely spaced simul- 
taneous observations. It should be noted that the assumption of independence 
of the gravity wave component profiles G would not be valid in this case, and 
this type of diffe^pnce yields a structure function which has significant var- 
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iation with the distance between the places of evaluation of the two profiles. 
However, since. the effects of tides and other phenomena are not removed, the 
horizontal structure functions yield valid results of truly irregular varia- 
tions only if the distance between profiles is kept small compared to the 
horizontal wave lengths of these tides and other phenomena. It should also 
be noted that the horizontal structure function used for velocity is actually 
the structure function of the horizontal velocity magnitude 

D(r) ■ {<[u(x + r) - u(x)] 2 > + <[v(x + r) - v(x)] 2 >}/2 (26) 

where u is the eastward velocity component, v is the northward component, x 
is the horizontal vector location of one velocity evaluation, r is the vector 
displacement between the locations of velocity evaluation, and r is the mag- 
nitude of the vector r. The definition (2&) is dependent on the magnitude r 
and not the direction of the displacement and no separate consideration of 

longitudinal and transverse structure functions are required. (A longitudinal 

2 

structure function would be of the form <[u(x + Ax) - u(x)] > and a transverse 

2 

structure function would be of the form <[u(y + Ay) - u(y)] >) . 

Harmonic Analysis . Data from three sites of Meteorological Rocket Net- 
work (MEN) data were collected and a harmonic analysis was performed on the 
data for periods of 5 days to 1000 days. For a particular site, only data 
obtained at the same time of day were accepted in order that the solar rad- 
iation error would be minimized, A function of the following form was fitted 
to the data by a least squares process: 

F(t) = A + A- sin ( 2 tt t /P + <j>) (27) 

o 1 
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where t is the time, is the mean value of the parameter, A ^ is the ampli- 
tude for the oscillation of period P, and $ is the phase angle. As the period 
P was varied from 5 days to 1000 days, a periodicity spectrum was calculated. 

Harmonic analysis was performed also on five sets of MRN data for tidal 
components. In the past, tidal analysis has been performed usually on data 
obtained from a high density of rocket launches over a basic period of about 
48 hours. However, there are frequenctly data from single rocket launches over 
a period of 5 to 10 days before or after the basic 48 hours. A preliminary 
error analysis indicates that the error in the harmonic analysis of the tidal 
components can be reduced significantly if input data includes these scattered 
data in addition to the data of the basic 48 hours. Solar radiation correct- 
ions were applied also to each data set. Then each data set was smoothed by 
the use of a polynomial smoothing function over five kilometer height inter- 
vals. Then, a function of the following form was fitted to the smoothed 
data by a least squares process: 

F(t) « A o + A i2 sin ( 2TTt / 12 + $ 12 ) + A 24 sin ( 27Tt / 2 ^ + $24^ ( 28 ) 

where time t is measured in hours, A q is the mean value of the parameter, A ^ 
and A^^ are the amplitudes of the 12 hour and 24 hour period tides, and 

$24 are phase angles for those tidal components. (The actual least squares 
fit was done in terms of sine and cosine terms equivalent to (28).) The re- 
sults of the harmonic analysis were subtracted from the solar radiation 
corrected data to find residual values. These residual values were evaluated 
at altitudes of one kilometer increments from a height of 45 km to 60 km. 

Thus, profiles of irregular variations, or gravity waves, were obtained. These 
gravity wave profiles were used to compute time structure functions. 
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3. MAGNITUDES OF GRAVITY WAVE AND 
PLANETARY WAVE VARIATION 

Resolution of the planetary wave component . Characteristics of the mag- 
nitude, probability distribution and structure of irregular winds at chemical 

release altitudes (- 90 - 130 km) were first studied by Woodrum and 
(4 5 ) 

Justus * . In this earlier analysis significant variations with height of 

the rms daily difference winds were observed. However, no variation with 
number of days separation between profile pairs was observed, and so, accord- 
ing to the concepts discussed in the previous section, it was considered tKat 
planetary waves do not propagate to the 90 - 130 km level with sufficient 
amplitude to be detectable. 

In the first part of this study ^ a large number of data in the 25 - 
200 km height range from various sources were collected for analysis by the 
daily difference method. These data included Meteorological Rocket Network, 
grenade, pitot tube, falling sphere, meteor winds, and satellite measurements 
as well as chemical release data. (For a complete list of data sources and 
references see the earlier report ^.) For the first phase of this study the 
daily difference analysis was performed by computing mean square daily differ- 
ences and averaging over all days of separation between 1 and 15. 

For the present phase of the study daily differences were computed sep- 
arately for each day of separation in order to determine if planetary wave 
contributions could be detected. The results are shown in Figure 1. This 
figure shows the mean square daily differences in the 45 - 65 km region, 

plotted versus number of days separation from 1 to 15, as a time structure 

2 

function, i.e. 2 a values versus time. 
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At the earlier times (1 to 72 hours) in Figure 1 actual time structure 
functions of the residual winds and thermodynamic variables from the tida 1 
analysis are plotted. The tidal analysis is discussed more fully in Section 5. 

The time structure data from the residual winds of the tidal analysis repre- 

sent an average over five sets of data made up of from 22 to 35 profiles in 
each set spread over a time span of 10 to 20 days with from 13 to 24 profiles 

occuring in a basic 48 hours for each set. A distinct peak at a time of 3 

hours (corresponding to a wave period of 6 hours) was found in the density, 
temperature, and pressure (not shown) structure functions. This indicates a 
preferred gravity wave period of 6 hours during the time from which data were 
available. No corresponding peak at 3 hours in the structure function of the 
residual winds was found. The time structure functions at one hour are def- 
initely smaller than those at subsequent hours and can be taken as an indi- 
cation of the upper limit e of the contribution from error and small scale 
turbulence. This rms upper limit e is related to the structure function at 
one hour by £ s (D(l)/2) (cf. equation 17). This relation results in the 
following estimate for the upper limits of error in the 45 - 65 km height 
range: 1% for pressure and density, 0.8% for temperature (~ 2°K) , and 3.5 to 
4.7 m/s for winds. These estimates are only slightly smaller than those com- 
puted by Avera and Miers v for this height region. Their estimates were 
2.6°K for temperature and 6.5 m/s for winds. 

The daily difference structure functions in Figure 1 are reasonably con- 
tinuous with the structure functions of the tidax residuals considering the 
smaller amount of data available for tidal analysis and the large amount of 
data in the daily differences. The daily differences in Figure 1 are aver- 
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aged over three MEN sites (Ascension Island, Cape Kennedy, and Fort Greely) 
and over altitudes from 45 to 65 km with data from 1964 through 1969. The 
daily difference magnitudes definitely increase with time separation and it 
appears that the daily differences at one day separation may be used as a 
reasonable estimate of the gravity wave contribution because of their general 
agreement with the large values of the tidal residual structure functions at 
time scales larger than 5 hours (i.e. periods > 10 hours). The daily differ- 
ence magnitudes are fairly uniform over the 7-15 days time separations and 
so an average over this interval can reasonably be taken as a measure of the 
combined gravity wave and planetary wave (up to 30 day period) contributions. 

Planetary wave magnitudes were extracted from the 7-15 day time separ- 
ation data for each of the three MEN sites mentioned above by the method des- 
cribed by equation (19), except that instead of one daily difference at large 
day number n, an average of the daily differences between 7 and 15 days was 
used. The results are presented in Figure 2. Primary contributions to the 
planetary wave magnitudes of Figure 2 would come from planetary waves in the 
period range 14 - 30 days. Data points represent averages over 10 km altit- 
ude sections centered about the point. The eastward planetary wav** wind com- 
ponent magnitude is larger than the northward component magnitude except at 
the high latitude site. Significant latitude variation in the 25 to 65 km 
altitude region was found. A steady increase in planetary wave magnitude 
with increasing latitude was observed for the pressure, density, temperature, 
and northward velocity component. The eastward component was reasonably sim- 
ilar at all latitudes. Table 1 gives the names and geographical locations of 
the MEN sites. Comparison of these results with those obtained by harmonic 
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Height, km Height 




TABLE 1. MEN SITE LOCATIONS 


Site Name 

Ascension Island, A.F.B. 
Fort Sherman, Canal Zone 
Barking Sands, Hawaii 
Cape Kennedy, Florida 
White Sands, New Mexico 
Point Mugu, California 
Fort Churchill, Canada 
Fort Greely, Alaska 


Latitude Longitude 


7° 

59' 

S 

14° 

25’ 

W 

9° 

20 ' 

N 

79° 

59' 

W 

22 ° 

02 ' 

N 

159° 

47’ 

W 

28° 

27' 

N 

80° 

32 ’ 

W 

32° 

23' 

N 

106° 

29’ 

W 

34° 

23' 

N 

119° 

07’ 

W 

58° 

44’ 

N 

93° 

49 * 

W 

64° 

00 ' 

N 

145° 

44* 

W 
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analysis are presented in Section 6. 

Above 65 km there were insufficient data for analysis of separate sites 
and so all 'data were combined for averaging. In the 65 - 85 km region and above 
it was found that the daily differences at 7 - 15 days time separation were not 
significantly larger than one day differences, which indicates no measurable 
propagation of planetary waves to this altitude range. However, as discussed 
more fully in the next section, there is evidence for propagation of selected 
wavelengths of planetary waves into the 65 - 85 km region. 

Gravity Wave Statistics . In the previous report ^ it was shown that 
probability distributions of the gravity wave fluctuations could be determined 
from the distribution of daily differences. In the light of the discovery of 
significant planetary wave contributions in the 25 - 65 km region, only daily 
differences over one day should be used to determine the gravity wave distri- 
butions. Details were given in the earlier report on the methods of determin- 
ing the distributions and the corrections applied. Figure 3 shows an example 
of the measured distribution of gravity wave perturbations determined from daily 
differences over one day. Values expected from a Gaussian distribution are 
shown for comparison. 

Some properties of importance in the study of the probability distribut- 
ions of a property x are the mean x = <x> and higher moments. The variance 

is given by = <r(x - x)S> the skewness is a = <(x - x)^>/0*^ and 

- 4 4 

the kurtosis is 3 = <(x - x) >/a . Table 2 shows the mean values of the 

gravity wave statistical parameters in the 50 to 80 km altitude region as 

determined from distributions of the single day daily differences. In Table 

2, as elsewhere u is the eastward velocity component and v is the northward 

component. In the 45 - 65 km range the kurtosis values 3 are all larger than 

the expected Gaussian distribution value of 3, whereas in the 65 - 85 km 
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region all of the g values except for the temperature g are less than 3, These 

revised values of kurtosis are similar to, though slightly different from, the 

previous results ^ in this height range. There are no changes from the higher 

altitude values previously reported. 

An analysis was performed to determine if the high g values in the 45 - 

65 km region could be due to intermittent gravity wave disturbance. Let us 

presume the following situation: A certain fraction y of the total observations 

consist only of Gaussian distributed error with standard deviation a^. The 

fraction y is known as the intermit tency. The remaining fraction (1 - y) of 

-k 

the measurements consist of Gaussian distributed gravity wave perturbations 

(9) 

with standard deviation a superimposed on the background error. Mosley 

S 

has shown that the kurtosis of the combined set of observations is given by 

(3 = 3(1 - yH 2 )/(l - y + yH) 2 (29) 

2 2 2 

where the ratio H is a ^ /( a e + a g ) • Equation (29) reduces to g = 3, the 

Gaussian value, in the limiting cases of no gravity waves (a = 0 or y = 1) or 

8 

no error (a^ = 0 or y = 0) . Figure 4 shows isocontours of g on the H - y 
plane. From the error estimates discussed above and the gravity wave magnitude 
estimates it was concluded that H - 0.2 for pressure and density, H - 0.3 - 0.4 
for winds and H - 0.5 for temperature in the 45 - 65 km region. From Figure 4 
it is possible that the observed g values for pressure and density (from Table 
2) could have been produced by an intermittency factor of 0.7 to 0.9. However, 
the remaining observed g values (i.e. for temperature and winds) cannot be 
explained by intermittency. 


This presumes several superimposed wave inodes interacting more-or-less 
randomly . 
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Figure 3. Probability distribution of gravity wave east- 
ward wind component in the 45-65 km height region. Shaded 
histogram gives observed distribution and solid dots show 
expected values from gaussian distribution of same stand- 



TABLE 2. Average gravity wave statistics 
in the 45 to 85 km region. N is the num- 
ber of independent values, a is the skew- 
ness and B is the kurtosis. 




45 - 65 km 


65 

- 85 km 



N 

a 

B 

N 

a 

B 

p 

1203 

+0.28 

6.56 

126 

+0.35 

2.18 

p 

1213 

+0.16 

5.68 

126 

+0.55 

2.56 

T 

717 

-0.54 

15.55 

126 

-0.23 

4.48 

u 

1934 

-0.07 

5.69 

59 

+0.07 

2.42 

V 

2033 

-0.08 

4.47 

59 

+0.02 

2.29 
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Figure 4. Values of kurtosis 3 at various values of 
the ratio H and the intermittency y as determined 
from equation (29). 


24 


Gravity Wave Magnitudes . The daily differences over one day separat- 
ion were used to compute revised gravity wave magnitude" in the 45 - 85 km 
region and new values in the 25 - 45 km region. Because planetary scale 
waves were not detectable in the 95 - 115 km region^^ the previous data, 
determined from averages over 1-15 day separation daily differences, were 
retained as the best estimates of gravity wave magnitudes above 85 km. 

Figure 5 shows the height variation of the measured average gravity 

wave magnitude. The data in Figure 5 were obtained from 1964 - 1969 MRN 

data one day differences from all eight sites listed in Table 1 plus all 

of the other falling sphere, grenade, pitot tube, chemical release, and 

satellite data collected. Meteor wind data from Garchy, France^^\ Durham, 

New Hampshire^*^ and Kazan, Russia^^ are shown separately in Figure 5. 

These meteor wind results are residual winds after explicit resolution of 

(13) 

the tidal winds by the method of Groves . The error bars in Figure 5 
were computed by the method outlined in Appendix A. 

The meteor wind results agree qualitatively although not in detail 
with the average results from the other data. However, it should be noted 
that the meteor results from Garchy and Durham were each determined from 
one set of measurements over a few days whereas the other results were 
obtained from several profile pairs from times spread over several years. 

The two gravity wave velocity components appear to be isotropic over 
the full altitude extent of Figure 5 with the possible exception of a 
slight predominance of the eastward component near 140 km. 

Linearized gravity wave theory with no dissipation or reflection 

2 2 

predicts that the gravity wave kinetic energy density o q <u + v > remains 
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Figure 5. Height variation of the magnitude of the gravity 
wave wind components. The eastward component is u and the 
northward component is v. 
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constant with height. However, Figure 6 shows that considerable dissipat- 
ion and/or reflection takes place, on the average, throughout, most of the 

L 

atmosphere. One exception is the 55 to 75 km height region, which corres- 
ponds to the decreasing temperature gradient of the atmosphere, a region 
[ where gravity wave reflection is not expected. At 90 km in Figure 6 the 
solid dot represents the MRN and other data and the open circle represents 
k an average value including the meteor wind data. Between 75 and 135 km the 
| average gravity wave kinetic energy density varies with height as 

2 2 

P q <u + v > * exp(-z/z Q ) (30) 

where the "scale height" value z q is 7.8 km, in good agreement with the 

value determined by Kochanski^"^ in the upper portion of this altitude 

range. Between 25 and 55 km the value for z q is 9.4 km. Extrapolation 
2 2 

of P Q <u + v > to 10 km indicates that jet stream level perturbations of 
about 3 m/s would be required to maintain the same rate of dissipation or 
reflective loss between 10 and 30 km as that observed between 30 and 50 km. 

If tropospheric propagation from the surface to 10 km is assumed to be at 
constant energy, analogous to the mesospheric propagation, then surface 
perturbations of about 1.8 m/s could maintain the dissipation or reflect- 
ive loss of the extrapolated curve. The change of slope above 135 km is 
due to a change in the density variation at this altitude while the gravity 
wave wind magnitude remains approximately constant with height. 

The height variation of the irregular pressure, density, and temper- 
ature is shown in Figure 7. The data are expressed as variation relative 
to the mean atmospheric parameters. Error bars in Figure 7 were also com- 
puted by the method outlined in Appendix A. 
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Height, km 



Figure 6. Height variation of the kinetic energy density 
of the gravity wave winds. Z Q values are energy fl scale 
height" values for equation (30) . The solid dots are the 
MRN and other data. The open circle at 90 km is an aver- 
age value including the meteor wind data. 
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RMS Value ,% 

Figure 7. Height variation of the magnitude of the relative pressure p’/p 0 
density p T /p 0 , and temperature T f /T due to gravity waves, and the theoretical 
p T /p 0 curve of Miller , equation (33), based on observed gravity wave winds 
from Figure 6. 



Figure 7 shows that the relative temperature variation magnitude between 
40 and 70 km is significantly smaller than for the relative pressure or den- 
sity, which are about equal. Between 80 and 100 km the relative pressure and 
temperature variation magnitudes become about equal and smaller than the 
relative density variations. Above 100 km the relative density variations 
become smaller than those for the relative temperature. 

At first one might think that the Boussinesq approximation should apply 
to the pressure, density, and temperature variations and hence that the rel- 
ative pressure variations should be much smaller than the relative density 
and temperature. However, as shown by Dutton and Fichtl^^, f G r motions of 
vertical scale not small compared to the scale height, the magnitudes of all 
three parameters are comparable. It will be shown later that the vertical 
scales are comparable with the scale height. Therefore, the relative values 
of the irregular thermodynamic parameters are related by 

p’/p o = p7p o + T'/T o (31) 

and in the mean square they would be related by 

<(p'/p o ) 2 > = <(p7p 0 ) 2 > + <t’/t o ) 2 > + 2<(p7p o )(T'/T o )> (32) 

where the last term is basically the correlation between the density and tem- 
perature variations, which in the case of a single component wave train would 
be dependent on the relative phase between the density and temperature wave 
variation. The changes discussed above in the relative magnitudes of the 
thermodynamic variables of Figure 7 therefore indicate a changing phase 
relationship with height. The changing phase relationship could be due to 
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nonlinear intermodal interactions or to selective modal absorption of a com- 
plex gravity wave field. In either case it would be difficult to jhstify the 
changing relative magnitudes with the assumption of a single mode gravity 
wave field. 

A comparison of the height variation of the irregular winds and irregular 
density, also shown in Figure 7 gives good agreement with a theory developed 
by Miller^^, which predicts relative density variations proportional to r'ms 
velocity variation by a factor which varies with height. The theoretical 
values of p'/p o are related to the mean square irregular winds by the formula 

|p'/p 0 l = {[(y - l)/c 2 + (l/gT Q ) dT Q /dz] <u 2 + v 2 >} 1/2 (33) 

This theory is applicable in the strict sense only to single mode gravity 
waves. The U.S. Standard Atmosphere of 1962 was used to evaluate the speed 
of sound c, the acceleration of gravity g, the temperature T q and its height 
gradient dT^/dz. The ratio of specific heats y was assumed to be constant 
at a value of 1.4. The agreement between the observed p f /P 0 variation and 
that predicted by the theoretical formula is remarkably good, considering 
the uncertainties introduced by the use of standard atmosphere values and 
the application to what is apparently a multi-mode gravity wave field. 

However, it does appear from the results shown in Figure 7 that the theo- 
retical formula corresponds more nearly with the magnitude of the irregular 
temperature T'/T than the density P i /p q - This discrepancy would have been 
improved considerably, however, if the density value at 110 km had been 
somewhat larger than the value actually observed. 

The three sites. Ascension Island, Cape Kenned^, and Fort Greely had 
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sufficient data for daily differences over one day separation so that sep- 
arate evaluations could be made for these sites. Figure 8 shows the altit- 
ude variation of the rms gravity wave magnitudes at these locations. There 
is no apparent system which describes the variation with latitude at all 
altitudes, although at many heights the magnitudes of winds and temperature 
fluctuations are smaller for the Fort Greely site. The large increase in 
wind magnitude between 50 and 60 km is compatible with the constant kinetic 
energy in this region, as discussed above. 

The large amplitude of daily differences at the White Sands site pre- 
viously reported^ ^ were found to be associated with the daily differences 
over large day separations and not with the single day daily differences. 
Hence, these large amplitudes of wind variation are to be associated with 
the planetary scale components rather than the gravity wave components. 
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Figure 8. Height and latitude variation of the gravity wave component. 
Ascension Island data is shown by the solid dots. Cape Kennedy by the squares, 
and Fort Greely by the triangles. 




4. SCALES OF GRAVITY WAVE VARIATION 

Vertical Scales . The vertical structure function of a parameter F(z) is 
given by 

D p (C) = <[F(z + 0 - F(z)] 2 > 

and, as shown in equation (23), the vertical structure function of the gravity 
wave components can be estimated from one half the vertical structure function 
of the daily differences AF^ taken over one day separation. Also, as shown by 
equation (25) the vertical structure function of the planetary wave components 
can be evaluated from the difference between the vertical structure functions 
of large day daily differences AF^ and single day daily differences AF^. 

Vertical scales of the gravity wave components can be estimated from the 
single day daily difference vertical structure functions by the technique 
discussed in Section 2. Figure 9 shows measured average vertical structure 
functions of gravity wave winds and temperature variations in the 45 - 65 km 
height range. Solid curves are fit through the data as determined from the 
exponential correlation functions suggested in Section 8. 

Figure 9 also shows the average vertical structure function of the large 
day daily differences. The difference in the two curves would be the struct- 
ure function of the planetary scale waves alone. It appears that there is 
some contribution to the structure functions from planetary scale waves of 
vertical wavelengths in the 10 to 30 km range, somewhat larger than planetary 
scale vertical wavelengths at lower altitudes . However, as the size of the 
error bars indicate, the firm conclusions which can be made on the vertical 
structure of planetary waves are somewhat limited. The discrete vertical 
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range determined from the single day daily differences (solid dots) and the vertical structure 
function of the 7-15 day daily differences (open circles) . Curves through the data points are 
computed from D(Az) = 2a 2 [1 ~ p(Az)] where the correlations are the exponential functions 
discussed in Section 8. 





wavelength reported earlier^ which appeared in data from some, but not all, 
sites seems to be associated with the large day daily differences, i.e. plane- 
tary scale waves, rather than the gravity wave components. 

The amount of data in the 65-85 km height range for vertical structure 
analysis was somewhat limited. Nevertheless some structure function results 
were obtained and these are shown in Figure 10. The vertical structure fun- 
ctions of the 10—12* day daily differences (i.e. up to 24 days period) are 
distinctly different from the one day differences and seem to indicate a 
preferred vertical wave length for the planetary waves of 20 km (structure 
function maximum at 10 km). However, the fact that the mean square daily 
differences were not significantly larger than the mean square single day 
differences, is taken to indicate a rather selective or sporadic propagation 
of these planetary waves into the 65-85 km region. 

Vertical structure functions above 85 km were not revised and remain 
unchanged from the results presented earlier^^. 

Horizontal Scales : There are no useable data available from which to 

determine horizontal structure of pressure, density, or temperature. The 
horizontal structure function of the winds in the 50 to 65 km region was 
slightly revised and is shown in Figure 11 (with a plotting error of the 
original graph corrected) . Also shown in Figure 11 is a plot of the hori- 
zontal structure function of winds in the 80-140 km height range. The solid 
curves fitted through the data in Figure 11 were evaluated from the engineer- 
ing approximation functions discussed in Section 8, using a values determined 
from the daily difference analysis. 
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6.4 m/s 
5 km 


Horizontal Winds 
50-65 km 


: 40 m/s 


200 km 


Horizontal Winds 
80-140 km 


Figure 11. Horizontal structure function of winds in the 50 to 140 
height region. Curves through the data points are determined from 
D(r) = 2o*[l - p(r)]with p(r) given by p(r) = [Rn( r l) + R22 (^i)]/ 2 
and the correlations R-j^ and R 22 given by equations (40) and (41) 
in Section 8. 







5. TIDAL HABMONIC ANALYSIS 


Tidal harmonic analyses have been performed previously by Justus and 
Woodrum^ on nine different sets of MRN data as an interim step to obtain- 
ing profiles of irregular winds, temperature, pressure and density. How- 
ever, the analyses were performed on the raw data which were not smoothed 
or corrected for solar radiation. Also, the length of the data records 
included only a basic 24-to-48 hours interval and did not take advantage 
of any additional data at scattered times around this basic interval. 

Thus, much better harmonic analysis results could be obtained if the above 
problems were corrected. Then, better profiles of the irregular components 
would be obtained, which, in turn, would give better estimates of gravity wave 
amplitudes by means of the time structure function. The revised time struct- 
ure function data were shown in Figure 1 for comparison with daily differ- 
ence analysis results. 

Smoothing and Solar Radiation Correction . First, the temperature data 

were corrected for solar radiation by a procedure given by Hoxit and 

( 18 ) 

Henry v . This procedure was valid only for certain types of temperature 
sensors. Hence, all data sets which used sensors other than these parti- 
cular types were eliminated from the analysis. It was found that five 
good data sets were left: one set from Ascension Island, A 4 F.B. around 
April 12, 1966; three sets from White Sands, New Mexico, around February 
7, 1964, November 22, 1964 and July 1, 1965; and one set from Fort Church- 
ill, Canada, around September 8, 1966. After the temperature data were 
corrected, the pressure and density data were recalculated from the hydro- 
static equilibrium equation and the equation of state. 
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Then, all the wind and thermodynamic data were smoothd by the use of a 
polynomial smoothing function over a height interval of five kilometers. 

Extended Time interval . Preliminary error analysis by Woodrum and 
(19) 

Justus indicated that significant differences in evaluated tidal paramet- 
ers may result if a data record of longer than 24 - 48 hours is used compared 
to results computed from a 24 - 48 hours data record. Further investigations 
have substantiated these earlier results. A preliminary analysis of errors 
was made by the following procedure: Thirty sets of fictitious data were 

created by superimposing a random error on smooth tidal data which included 
the periods of 24 and 12 hours as given by the following equation. 

Y (t) = 10 sin + 5 sin + 6) + 12R (34) 

where R is a random number between -1 and 1. The amplitude of the random com- 
ponent in equation (34) is representative of the fluctuations produced by 
gravity waves and traveling planetary scale waves of a few days period, as 
determined in Section 3. 

Harmonic analyses on the 30 data sets were performed and standard deviat- 
ions of the amplitudes and phases were found. Then, the time length of each 
data set was reduced and the harmonic analysis was repeated. This process 
gives the standard deviation as a function of the time length of the data sets 
as shown in Figure 12 and 13. These figures show results of two analyses for 
diurnal and semidiurnal oscillations where the data sets first had values at 
every hour and then had values only for every third hour. It is seen from 
the figures that the errors seem to be rather small unless the time length 
of the data set is less than about three days. The basic time length of the 
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Figure 12. Amplitude Error Analysis of Simulated Data. 
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Figure 13. Phase Error Analysis of Simulated Data. 
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observed MRN data sets used in this report are of the order of two days. Thus, 
Figure 12 and 13 show that the harmonic analysis of these data should be rea- 
sonably accurate but that the time length of the data is short enough to be 
questionable. However, the data values used in these experimental data are 
not evenly spaced in time as is the case for the ficticious data. In order 
to better approximate the actual data, ficticious data values were selected 
at times which closely corresponded to the times of the actual data. Harmonic 
analyses were performed on these ficticious data and standard deviations were 
calculated as shown in Table 3. The smaller number of data points for each 
site corresponds to the basic time interval whereas the larger number of 
data points corresponds to the extended data. 

In comparing the standard deviations in Table 3 of the basic and extended 
time intervals, one sees that the standard deviations may be reduced in some 
cases as much as 25% to 30%. The magnitude of these reductions are very sig- 
nificant . 

An example of the comparison of an harmonic analysis of the experimental 
data for the two time intervals is shown in Figure 14. This figure shows the 
northward wind component from White Sands, set 1. The basic time length of 
the data was 2 days. Although the two data sets give similar results, there 
are definite differences. The amplitudes of the basic time interval are 
about 25% smaller at the higher heights and the phase is much more erratic 
than that of the extended time interval. These results indicate that one 
should definitely include the additional data scattered close in time 
around the basic time interval in the harmonic analyses of the atmospheric 
data. One must be careful, however, to keep the extended time interval short 
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TABLE 3. COMPARISON OF STANDARD 
DEVIATIONS OF BASIC AND EXTENDED DATA 


Number of Diurnal Diurnal Semidiurnal 

Data Points Wind Wind Phase Wind 

(m/sec) (Deg.) (m/sec) 


Ascension Island 


24 

1.38 

10.92 



2.55 

32 

1.33 

10.19 



2.47 



White Sands, 

Set 

_1 


20 

2.51 

13.26 



2.12 

35 

2.03 

10.33 



1.92 



White Sands, 

Set 

_2 


17 

2.12 

18.76 



2.62 

28 

1.55 

12.01 



2.08 



White Sands, 

Set 

_3 


13 

2.73 

18.03 



2.29 

25 

2.48 

15.45 



1.98 



Fort Churchill 



19 

1.93 

13.76 



2.16 

22 

1.99 

12.49 



2.15 


Semidiurnal 
Wind Phase 
(Deg.) 


57.67 

47.07 


30.64 

22.39 


45.27 

37.11 


36.01 

24.92 


48.66 

32.75 
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Figure 14. Comparison of Harmonic Analysis of Basic and Extended Data 
for the Diurnal Northward Wind at White Sands, Data Set 1. 
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enough so that seasonal changes in tidal phase or amplitude would be small. 

Time Structure Analysis . Harmonic analyses were performed on the five 
sets of data with the above mentioned corrections. The results were sub- 
tracted from the data which were corrected for solar radiation only; and 
thus, irregular data profiles were obtained. Then these irregular profiles 
were used in the time structure function to give an estimate of the gravity 
wave amplitudes as shown in Figure 1. The harmonic analysis results compared 
favorably with previous work that had been performed on essentially the same 
d at a (20 - 25 >. 

Compatibility of Wind and Thermodynamic Data . A major difficulty in 
performing the harmonic analysis of upper atmospheric data is the insuffic- 
ient quantity of properly space time series data. However, the need for the 
large quantity of data could possibly be reduced if the number of independent 
variables in the tidal equations of motion could be reduced. 

In theory it is possible to infer temperature, pressure, and density 
data from the wind data^^’^^. In general, the process involves substitut- 
ion of the wind component data into the east-west equation of motion to ob- 
tain geopotential height variations at several pressure levels. Then, these 
height variations are converted into pressure, density, and temperature var- 
iations. This method would make the pressure, density, and temperature var- 
iations dependent on the wind variations and thus, would certainly reduce 
the number of independent variables in an harmonic analysis. However, the 
method is very sensitive to errors in the wind components. To check the 
feasibility of applying this method to the data in the present report, we 
used the method to calculate the temperature, pressure and density variations 
from the observed wind variations for one data set from White Sands and then. 
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compared these results with the observed thermodynamic data. The comparison 
of the two results was very poor. Similar results were obtained by Beyers et 
al.^^\ Consequently, the conclusion of this analysis is that this method 
should not be used as a means of reducing the number of independent variables 
because it seems to be too sensitive to errors. 

Another means to reduce the number of independent variables is to restrict 
the phase difference between the northward and eastward wind components to 
90°. Theory predicts that if the tidal motion contains only one mode of 
oscillation, the zonal and meridional wind components will be in phase quad- 
rature. However, if more than one mode of oscillation is present, the phase 
relationship between the zonal and meridional wind could take on any value. 
Inspection of the results of the tidal analysis from the five data sets in 
the present report indicates that there is no consistent phase relationship 
between the zonal and meridional wind components. Thus, this indicates that 
more than one mode is present in the tidal motion and the phase difference 
restriction to a value of 90° would be invalid. 

Thus, the conclusion is that, at present, there is no means to by-pass 
the need for the large quantity of properly spaced data in order to perform 
the tidal harmonic analysis. 
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6. LONG PERIOD VARIATIONS 


Sufficient MRN data were found for three sites to perform harmonic analy- 
sis using periods of 5 days and longer. The data were chosen for each site 
such that the data values occurred at the same time of day, plus or minus ten 
minutes, for a time interval covering about six years, 1964 through 1969. The 
sites and the local times of day were Ascension Island, A.F.B., 1345 hours; 
Cape Kennedy, Florida, 1000 hours; a n d Fort Greely, Alaska, 1000 hours. Site 
latitudes and longitudes are given in Table 1. 

No solar radiation corrections were made on the temperature data since 
the corrections for solar angle variation at the same time of day throughout 
a year are quite small. The radiation errors would be essentially uniform in 
the data and would not affect the harmonic analysis. The only non-uniform 
radiation errors in the temperature data would be errors due to seasonal var- 
iations of the solar angle which are relatively small (less than 2°K or about 
0 . 8 %). 


The number of data points for the analysis vary with height from a maximum 
of 409 to a minimum of 65. The distribution of these data with seasons of the 
years are as follows: 




Winter 

Spring 

Summer 

Fall 


Ascension Island 

26% 

18% 

30% 

26% 


Cape Kennedy 

30% 

22% 

22% 

26% 


Fort Greely 

24% 

35% 

25% 

16% 

The 

data are rather evenly 

distributed by seasons. 




Harmonic analysis was 

performed by 

fitting one 

period at a time 

to the data 

and 

thus, a periodogram was obtained. 

The heights 

included in the 

analysis were 


45 km through 60 km in increments of 5 km. For example Figures 15 through 23 


*8 
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Figure 15. Periodogram of Wind Components for Ascension Island 
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Figure 16. Periodogram of Wind Components for Cape Kennedy. 



Figure 17. Periodogram of Wind Components for Fort Greely. 
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Figure 21. Periodogram of Pressure and Density for Ascension Island. 



Figure 22. Periodogram of Pressure and Density for Cape Kennedy. 



Figure 23. Periodogram of Pressure and Density for Fort Greely 




show the periodiograms for the winds, temperature, pressure and density at the 
height 45 km. The predominant oscillations shown by these figures are the 
annual, semi-annual, and quasi-biennial (period = 780 days) oscillations. The 
temperature seems to have a strong oscillation at the 600 days period also. 
Notice that the long period northward wind is very small compared to the east- 
ward wind both at Ascension Island and Cape Kennedy and only becomes signifi- 
cant at the higher latitude of Fort Greely. 

Quasi-biennial Oscillations . The quasi-biennial oscillation appears as 
definite peaks on. the periodograms of the winds, pressure and density. The 
temperature curves do not show any preference for the quasi-biennial period 
as compared to adjacent periods. Table 4 shows the amplitude variations of 
this oscillation with height and latitude. Generally, the oscillation de- 
creases in amplitude with increasing heights near the equator and increases 

im amplitude with increasing height at the higher latitudes. These results 

(27-35) 

are in agreement with previous work . Figure 24 shows the month of 

maximum amplitude of the quasi-biennial oscillation for the different sites. 
The abscissa of this figure has two scales. The top scale shows the month 
of maximum amplitude in the years 1966-1968 as calculated. The bottom scale 
shows the month of maximum amplitude for the years 1972-74 assuming that the 
oscillation remains coherent. 

Annual Oscillation . The annual oscillation is definitely the most pre- 
dominant of all the planetary waves as shown by the periodograms. Figure 
25 shows the month of maximum amplitude for the different sites. Table 5 
shows the variations of the amplitudes with height and latitude. The wind 
amplitudes for each site are rather constant with height. The eastward wind 
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Height, km 
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IVlonth of Maximum 

Figure 24. Month of Maximum for Quasi-Biennial Oscillations. 
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Month of Maximum 


Figure 25. Month of Maximum for Annual Oscillations. 
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TABLE 4. AMPLITUDE VARIATION WITH HEIGHT 
AND LATITUDE OF QUASI-BIENNIAL OSCILLATION 


Height 

(km) 

Ascension 

Island 

Cape 

Kennedy 

Fort 

Greely 


Eastward Wind 

(m/sec) 


45 

19.2 

11.9 

5.9 

50 

15.5 

15.9 

7.6 

55 

5.9 

16.2 

11.4 

60 

1.5 

34.5 

16.4 


Northward Wind 

(m/sec) 


45 

0.55 

1.15 

3.92 

50 

1.17 

0.78 

3.86 

55 

1.19 

0.35 

7.21 

60 

2.62 

4.65 

14.87 


Temperature 

(K deg) 


45 

1.86 

5.97 

3.56 

50 

1.87 

4.46 

4.08 

55 

0.64 

8.03 

4.13 

60 

0.84 

17.89 

6.45 


Pressure (% of mean) 


45 

1.5 

0.3 

5.6 

50 

1.4 

0.8 

7.0 

55 

1.2 

1.2 

7.6 

60 

0.9 

1.9 

8.4 


Density (% o 

f mean) 


45 

1.5 

0.8 

5.0 

50 

1.5 

0.5 

6.5 

55 

1.3 

0.9 

7.5 

60 

1.8 

1.5 

6.5 


55 





TABLE 5. AMPLITUDE VARIATION OF ANNUAL 
OSCILLATION WITH HEIGHT AND LATITUDE 


Height Ascension Cape Fort 

(km) Island Kennedy Greely 


Eastward Wind (m/sec) 


45 

18.4 

37.0 

23.8 

50 

16.6 

42.6 

26.8 

55 

17.8 

43.0 

29.0 

60 

14.4 

40.1 

22.3 


Northward Wind (m/sec) 


45 

2.9 

3.2 

13.9 

50 

2.9 

2.8 

15.8 

55 

3.1 

2.2 

15.1 

60 

3.0 

2.3 

14.6 



Temperature (K deg) 


45 

0.7 

7.0 

21.1 

50 

1.1 

4.3 

16.6 

55 

1.8 

7.8 

12.6 

60 

2.9 

21.7 

9.6 



Pressure (% of mean) 


45 

1.4 

5.1 

20.9 

50 

1.0 

4.9 

25.5 

55 

0.5 

4.8 

28.7 

60 

0.7 

3.6 

28.4 


Density (% of mean) 


45 

1.9 

5.9 

12.8 

50 

1.2 

5.1 

19.2 

55 

1.1 

4.7 

23.7 

60 

1.7 

4.0 

24.7 
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has its largest oscillations at Cape Kennedy. The northward wind is insigni- 
ficant for the two lower latitude sites but becomes appreciable at Fort Greely. 

The temperature oscillations, which are in units of K deg, and the pressure 
and density oscillation, which are in units of per cent of the mean, behave in 
similar fashions. Their amplitudes are small near the equator (less than 3) 

but increasing with latitude to an amplitude in the 20* s at Fort Greely. These 

(34) 

oscillations agree reasonably well with the results of Cole 

Semiannual Oscillation . The semiannual oscillations seem to be subordi- 
nate to the annual oscillations everywhere except near the equator. At the 
site of Ascension Island (especially in the eastward wind) , the semiannual 
amplitudes are larger than the annual amplitudes at the lower heights. How- 
ever, the semiannual amplitudes decrease rapidly with height so that the 
annual oscillations appear to dominate at the higher heights even near the 
equator. 

Figure 26 shows the month of maximum amplitude for the different sites. 
Table 6 shows the amplitude variations with height and latitude. The eastward 
wind, as stated above, seems to decrease im amplitude with increasing height, 
but shows a maximum at the height of 50 km which agrees with experimental re- 
suits of Reed'' . The eastward wind amplitude remains essentially constant 
with height at Cape Kennedy and seems to increase with height at Fort Greely. 

The northward wind amplitude again is insignificant for the two lower latitudes 
but does become significant at Fort Greely. 

Again, the temperature, pressure and density oscillations act similarly 
by having small amplitudes near the equator but increasing with latitudes to 

/o i 

amplitudes of 10 to 15 at Fort Greely' ’ 
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Height, km 



Month of Maximum 

Figure 26. Month of Maximum for Semiannual Oscillations. 
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TABLE 6. AMPLITUDE VARIATION OF SEMIANNUAL 
OSCILLATION WITH HEIGHT AND LATITUDE 


Height Ascension Cape Fort 

(km) Island Kennedy Greely 


Eastward Wind (m/sec) 


45 

24.9 

15.7 

8.3 

50 

25.8 

17.9 

11.7 

55 

18.5 

17.7 

15.2 

60 

8.9 

18.2 

19.9 


Northward Wind (m/sec) 


45 

0.3 

2.4 

5.9 

50 

0.5 

2.9 

5.2 

55 

1.0 

2.9 

4.3 

60 

1.9 

2.2 

2.6 



Temperature (K deg) 


45 

2.2 

6.8 

8.3 

50 

1.8 

3.8 

7.9 

55 

1.4 

7.0 

7.2 

60 

2.7 

22.8 

4.8 



Pressure (% of mean) 


45 

1.7 

1.5 

8.4 

50 

1.8 

1.4 

10.5 

55 

1.4 

1.3 

13.7 

60 

1.1 

0.7 

15.8 



Density (% of mean) 


45 

1.3 

2.4 

6.3 

50 

1.5 

1.8 

8.0 

55 

1.8 

1.0 

10.9 

60 

1.0 

0.8 

14.3 
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Shorter Period Oscillations. Theories have been presented ^ which 

predict that planetary waves of periods 5 to 50 days may penetrate upward into 

the stratospheric regions. A few observations have been made which seem to 

(41 42) 

substantiate these theories * and recent refinements to the theory have 

(43) 

been offered . In the present analysis, attempts to identify planetary waves 
of periods 5 to 50 days from the periodograms were made. There seem to be indi- 
cations in the periodograms of oscillations with these periods; however, the 
oscillations are not consistent for the different sites or for the different 
parameters. For example, oscillations with periods of about 6 days and 25 days 
appear to peak in the eastward wind at Cape Kennedy with an amplitude of about 
10 m/sec, but not at the other two sites. In addition, the other parameters 
do not show peaks at these two periods except possibly the pressure and density 
oscillations at Fort Greeley. 

The periodograms indicate that the largest oscillations in the wind and 
temperature in the periodic range of 5 to 50 days occur at Cape Kennedy with 
amplitudes of about 10 m/sec and 7 K deg. The other two sites show very small 
oscillations in the winds and temperatures for these periods. The oscillations 
in pressure and density seem to be largest at Fort Greely with amplitudes of 
about two per cent of the mean and seem to be smallest at Cape Kennedy. 

The amplitudes of the 14 to 30 day periods were analysed also by the meth- 
od of daily differencing (See Section 3) . These amplitudes are shown for com- 
parison in Figures 15 through 23 by the arrows on the left side of the figures. 
In every case except two, the amplitudes as shown by the periodograms are 
much smaller than those predicted by the daily difference method. The prob- 
able cause for the discrepancy is that the harmonic analysis method assumes 
that the oscillations remain coherent for the entire time of investi- 
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gat ion. It is highly unlikely that a planetary wave of period 5 to 50 days 
would remain coherent for six years. Thus, if the waves were incoherent, the 
harmonic analysis method would underestimate the amplitudes of these waves. 

To do a harmonic analysis for these planetary waves, one needs to divide the 
data into smaller time intervals so that coherence could be attained. Then, 
the results of all the time intervals could be averaged. 
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7 . SOLAR CORRELATIONS 


Correlation at Chemical Release Heights. A significant correlation has 
been found between the northward wind component and the F^q ^ solar flux 
(or Zurich sunspot number R^) in the altitude range 110 - 140 km. Density 
and temperature variations near the 100 km level have also been correlated 
with solar activity The eastward wind component, however, was not found 
to be correlated with solar activity indices in a significant way. The north- 
ward wind correlation was such as to give a regression relation 

V N - 42 + 0.85 F ltK7 (35) 

Since the average value of F^ ^ is approximately 80, a 10% change in F^q ^ 
would cause a change of roughly 7% in the northward wind component. 

It is believed that the dependence of northward velocity on F^q ^ expressed 
by equation (35) is the result of the generation of gravity waves by magnetic 
storms and their subsequent propagation equatorward. The southward propagating 
gravity waves would produce a predominately north-south velocity perturbation. 
The generation of gravity waves by magnetic substorms was proposed earlier by 
Testud and Vasseur*^^. 

Observations of increased eastward wind components at higher altitudes 

(47) 

(near 160 km) during a geomagnetic storm were reported by Smith but the 
increased eastward component was suggested as being due to a thermal wind 
relationship with increased polar temperature over equatorial temperature. 

A solar cycle dependence of prevailing winds (both northward and eastward 
components) and semidiurnal component winds have been reported by Springer and 
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Schminder^ 8 ^ . These observations at about 95 km altitude indicate an in- 
crease in the semidiurnal component with solar activity, but a decrease of 
the prevailing winds with solar activity. 

It was decided that, rather than a correlation analysis, a superposed 

(49) 

epoch analysis, similar to that applied by Woodbridge at the 300 mb level, 

would be most appropriate. The superposed epoch analysis would consist of 
mean square values of a parameter (e.g. gravity wave or planetary wave wind 
or thermodynamic variables) plotted versus number of days before or after a 
key day, i.e. a day of high solar activity. However, it was considered that 
at this time the amount of data above MEN height, is such that too few data 
would be from times close to key days for such an analysis to be performed 
(See below for discussion of the method of key day selection) . Thus the super- 
posed epoch analysis above 65 km was postponed until more data can be added to 

t 


that already accumulated. In the meantime it is felt that the correlation given 
by equation (35) can be used to estimate the gravity wave magnitude through 

AV N * °' 85 A ^10. 7 * 

Correlation at MEN Heights . The superposed epoch analysis was employed on 

the MRN data from three sites: Ascension Island, Cape Kennedy, and Fort Greely. 

The first problem was to define the key days to be used in the analysis. Noon- 

kester^"*^ used as a definition any day in which the international character 

(49) 

figure C^ increased by 1.0 or more from the preceding day. Woodbridge con- 
sidered as a key day any day on which the planetary geomagnetic index was 
15 or larger and had first reached this value through a one-day rise AA^ 
greater than 11. The Noonkester definition yielded 56 key days in the 1964-69 
period studied, or about one every 39 days on the average. The Woodbridge de- 
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finition yielded 151 key days or about one every 14 or 15 days on the average. 
Therefore a modified definition was developed to reconcile these differences. 
The following definition was adopted: A key day is one in which is greater 

than 15 and first became greater than 15 through either (1) a one-day change 
AA > 11 and on which also the one-day change AC. > 1, or (2) a one day change 

p - l - 

AA > 18 regardless of the AC. value. With this definition 50 of the 3b 
P - 1 

AC^ > 1 days corresponded to AA^ > 11 days and there were 28 additional 
AAp > 18 days which did not correspond to AC^ > 1, with a total of 78 key days. 
Table 7 gives a yearly and monthly breakdown of the occurrence of key days 
determined by this method. 

After the key days had been selected, daily difference magnitudes in the 
25 to 45 km and 45 to 65 km height ranges were averaged with respect to number 
of days before or after a key day (based on the date of the earlier profile of 
the pair being differenced). Daily differences over one day separation, to 
check for gravity wave dependence on solar activity, and over 7 to 15 days 
separation, to check for planetary wave dependence on solar activity, were 
tested. In all cases, in both height ranges, no significant variation of mag- 
nitude with respect to days after a key day was determined. 

Thus it appears, on the basis of this study, that the solar correlation 

(44 45) 

found near the 100 km level * ' and that found at lower stratospheric 

(49-52) 

heights are not connected by processes which propagate through the 25 - 

65 km region. This supports the idea that the correlations near the 100 km 
level are due to direct geomagnetic storm generation of gravity waves , 
while the lower stratospheric, and possibly surface correlations are due to 
a complicated chain of x-ray nucleation, cloud formation, and surface rad- 
iation balance effects, such as proposed by Roberts . 
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TABLE 7. YEARLY AND MONTHLY DISTRIBUTIONS 
OF THE 78 KEY DAYS 


Year 

Number 

Month 

Number 

Month 

Number 

1964 

17 

Jan. 

6 

July 

6 

1965 

7 

Feb. 

12 

Aug. 

5 

1966 

15 

Mar. 

9 

Sept . 

9 

1967 

16 

Apr. 

6 

Oct. 

6 

1968 

11 

May 

9 

Nov. 

2 

1969 

12 

June 

4 

Dec. 

4 
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It will be possible, with the addition of a few more data sets, to do a 
superposed epoch analysis of winds in the chemical release altitude region. 
With this approach it should be possible to compute gravity wave magnitudes 
during geomagnetic storms and during non-storm conditions. 
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8. PROCEDURE FOR EVALATION OF EFFECTS ON 
SPACECRAFT DUE TO ATMOSPHERIC VARIATION 

The analysis of structure functions indicate that the gravity wave per- 
turbations in wind, pressure, density, and temperature in the 25 - 200 km reg- 
ion have a continuous spectrum. Thus the effects of gravity waves on a space- 
craft traversing these altitudes can be analyzed in terms of the spectral 
techniques of turbulence. 

Planetary waves in the 45 - 65 km region exhibit single modes of specific 
vertical wavelength (above 37 km) and in the 65 - 85 km region a preferred 
vertical wavelength of about 20 km is found for the planetary scale waves. 
Atmospheric tides also exhibit discrete wavelengths. Statistical but non- 
spectral methods for considering these discrete wavelength phenomena are pre- 
sented later in this section. 

The coordinates of the analysis are indicated in Figure 27. The x, y, 
z system has x and y axes in the horizontal plane and z in the vertical dir- 
ection. The x direction has been conveniently chosen to be along the azimuth 
of the spacecraft trajectory. The elevation of the spacecraft trajectory is 
the angle © and the x’, y 1 , z f axes are those which are rotated through this 
angle 0. Thus x 1 is along the direction of flight of the spacecraft. The 
displacement r along the x f axis is the distance of travel of the spacecraft 
in an arbitrary time interval t. With the speed of the spacecraft given by 
U, then r and t are related by r = Ut. The displacement r can be broken into 
component displacements r^ and r^ in the horizontal and vertical directions 
given by 

r x - U cos 6 t r 3 = U sin 0 x (36) 
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displacement r is along trajectory path. 


Correlation Functions. Tne structure functions measured and reported 


both in section 3 and 4 and in the earlier report ^ can be related to corre- 
lation functions, which in turn can be related to spectra. The measurements 
indicate that both scalar and vector spectra should have a high frequency 
range power law behavior with exponents ranging from -1.6 to -2.5. The high 
frequency range exponent for the von Karman spectrum is -1.67 and for the 
Dryden spectrum it is -2.0. The Dryden spectrum, which is based on a linear 
exponential correlation function, is often used in place of the von Karman 
spectrum for atmospheric turbulence and the differences between the two have 
little engineering significance because the major differences occur only at 
the relatively unimportant high frequencies. Since the Dryden exponent falls 
roughly in the middle of the observed range of exponents then the Dryden spect- 
rum can be used in all cases without serious engineering error. However, the 
anisotropy of the irregular variation fields must be taken into account by the 
use of different vertical and horizontal scales. 

For engineering purposes the following spatial correlation functions can 
be used in the 25 to 200 km altitude range: For scalar quantities (pressure, 

density, and temperature) use 

R(r) = e~ r/L (37) 

where r/L is given by 

(r/L) 2 = (r 1 /L 1 ) 2 + (r^) 2 (38) 

and and are respectively the horizontal and vertical scales given in 
Table 8 for gravity wave pressure, density, and temperature variations. Sub- 
stitution of (36) into (38) provides a relation between the length scales 
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and the trajectory angle 0 


L 


L^L^/ (L^ 2 sin 2 


, t ’2 2 

+ cos 


e) 


1/2 


(39) 


For spatial correlations of velocity components use 


*u ( ? 

- (1 - 

r 3 2 /2rL u ) e r/Lu 

(40) 

R 22 ( r ) 

- (1 - 

r/2L ) e“ r/Lu 
u 

(41) 

R 33 ( ! } 

- (1 - 

r i 2/2r V e ' r/L " 

(42) 

R 13 ( ! ) 

+ R 31 (r) 

■ < r l r 3 /r V e_r/L “ 

(43) 


where Rj_(r) is the correlation <u^(x) u^ (x + r)>/a^Oj, is 
the ith velocity component, and and are defined by 


the rms value of 


(r/L u ) 2 - (t./L^) 2 + (r 3 /L u3 ) 2 (44) 

2 2 2 2 1/2 

L = L L 0 /(L / sin 0 + L 0 cos 6) 1 (45) 

u ul u3 ul u3 

(r/L w ) 2 - (rj/L^) 2 + (r 3 /L w3 ) 2 (46) 

L „ ‘ L ul L „3 /a „l 2 s1 " 2 9 + h J c ° s2 9)1/2 (47) 


where L L 0 , L ,, and L ^ are horizontal and vertical scales given in Table 
ul* u3* wl w3 

8, along with the appropriate mean a values for the irregular horizontal winds 
u and v (u^ and u^) and vertical component w (u^) . 

Scalar and Velocity Spectra . Spectra computed from these correlation 




functions are as follows: For the spectrum of scalar quantities use 


2a 2 L/ir[l + (kL) 2 ] 


(48) 


4>(k) * 

where a is the rms value of the quantity whose spectrum is being evaluated. 

The wave number spectra in (48) is for k in radians /km, with L determined from 
(39), or in terms of frequency u) in radians/sec use 

<p (cu) - 2a 2 L/ttU[ 1 + (u L/U) 2 ] (49) 

Table 8 gives values of a for gravity wave pressure, density and temperature 
variations. For trajectory angles 0 less than about 5° the spectrum computed 
by either (48) or (49) can be considered the same as for level flight, i.e. 
use for L. The trajectory angle 9 is less than 5° throughout the entire 
descent phase and above 85 km on the ascent phase of the typical space shuttle 
trajectory (see Table 9 for the typical values used here) . For the transverse 
velocity spectrum in the horizontal direction use 

* 22 '(k) = a u 2 L u [l + 3(k L u ) 2 ]/tt[ 1 + (k O 2 ] 2 (50) 

for the wave number spectrum with determined from (45) or in. terms of fre- 
quency use 

<f> '(u>) = a 2 L [1 + 3(w L /U) 2 ] / tf U[1 + (u L /U) 2 ] 2 (51) 

22 u u u u 

For trajectory angles less than 2° horizontal flight may be assumed and 
used for L u in (50 or (51) . For the transverse velocity spectrum in the ver- 
tical plane use 


71 


, a 2 cos 2 9 L [1 + sin 2 0 + (2 + cos 2 9)(k L ) 2 ] 

*33 00 - — W 


ir [1 + (k L ) 2 ] 2 
w 


a 2 sin 2 0 L [1 + cos 2 0 + (2 + sin 2 0) (k L ). 2 ] 
u u_ u 

it [1 + (k L u ) 2 ] 

2a a sin 2 0 cog 2 0 L [1 - (k L ) 2 ] 
u w w w 

TT [1 + (k L ) 2 ] 
w 


(52) 


For 6 < 2° horizontal flight may be assumed and (52) reduces to 


+ 33' 00 " °w 2 L »1 U + 3(k L wl )2]/,tl + (lt L »1 )2]2 (53) 


w r 


or, expressed as a frequency spectrum 

* 33 ' (a)) = a w 2 L wl [l + 3(a) L w 1 /U) 2 ]/ttU[ 1 + (a) L wl /U ) 2 ] 2 (54) 

Because of small values of a and L the transverse vertical spectrum com- 

w w 

puted by (52) remains a factor of ten or more smaller than the horizontal trans- 

t Q 

verse spectrum *22 even U P to trajectory angles of 15 . Therefore, the simple 
forms (53) or (54) may be used in all cases in which the vertical spectrum, be- 
cause of its smaller magnitude, would not have to be evaluated with extreme 
accuracy. 

The longitudinal spectrum would have a complicated angular form analogous 
to (52) but since it is no more angular sensitive that the transverse spectra, 
and since the longitudinal response modes are usually unimportant, only the 
level flight spectral form need be considered. With 0=0, the longitudinal 
spectrum is 

* u '(k) - 2o u 2 L u] Ml + (k L ul ) 2 ] (55) 
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Table 8. Horizontal and vertical scales and mean standard deviation values for irregular winds, pressure, 
density, and temperature. Values in parentheses indicate interpolations, extrapolations, or 
values assumed by analogy with other variables. 
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or in frequency form 


<t> 11 (u) 


2 a u 2 L /irU[l + (to L ul /U) 2 ] 


(56) 


TABLE 9. TYPICAL 

SPEEDS U AND 
Values were i 

TRAJECTORY 

determined 

ANGLES 6 FOR 
from Smith et 

SPACE SHUTTLE TRAJECTORIES. 
al< 53 \ 

Height 

Ascent 



Descent 

Range 





(km) 

u 

6 

U 

e 


(km/ sec) 

(deg.) 

(km/sec) 

(deg.) 

45 - 65 

2.4 

14° 

4.6 

0.6° 

65 - 85 

3.4 

7° 

4.7 

0.2° 

85 - 110 

5.9 

2° 

4.6 

0.8° 


Probabilities and Exceedance Values . No measurements have been made of the 
probability distribution of the standard deviation p(a). Only the average val- 
ues have been determined. The form of p( a ) which should be assumed is the same 
as used for turbulence at lower altitudes. 

p(a) = ( / 2 /tt b) exp (- a^/2b 2 ) (57) 

from which it can be shown that the average value of o is a - b / 2/tt, so that 
the values of b (the standard deviation of a) to be used are given by 

b - a/77 2 = 1.253 a (58) 

where a values can be obtained from Table 8. Therefore the number M of fluct- 
uations of a load quantity y that exceed the value y* is given by 
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M(y*) 


(59) 


= N q exp(- |y*|/1.253 a A) 


where N q is given by 


N 


(1/2* o y ) [ j” u> 2 ♦y(M) du>] 1/2 


( 60 ) 


in which <J> y is the spectrum of y and a y is given by 


r 


a « [ (j) (a)) d ( d] 


1/2 


(61) 


For a simgle element of the spacecraft the relation between the spectrum <p (to) 
of the response y and the spectrum <f> (to) causing the response is 


♦ (<*>) * ♦ (u)Ih (to) | 


(62) 


where H is the transfer (frequency response) function. The parameter A in 
equation (59) is given by 


A « Oy/o (63) 

(54) 

Thus a standard type of exceedance model such as that employed by NASA can 
be used with b, the standard deviation of a, given by equation (53). With fur- 
ther study perhaps p(a) could be divided into quiet or non-storm conditions and 
active or storm conditions based on mid-winter warming conditions, especially 
at the lower altitudes, and based on solar activity, especially at the higher 
levels. Perhaps a better relation than (57) could be determined by actual 
measurement of the probability of occurrence of a values. Seasonal variation 
of p(a) may also exist. 

Latitude and Other Variations. If it is desired to take into account the 
different irregular variation magnitudes at various latitudes then the average 
standard deviation values from Figure 8 may be used in the 25 - 65 km altitude 


range. 
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If account is to be taken of the relative amount of time spent in the var- 
ious flight regimes (latitude, altitude, etc.) then consider the case of a total 
of k mission segments, with time t^ spend in each mission segment, and a total 
flight time of T. The total exceedance rate for all mission segments is 


M(y*) = 

and the paraMters v 

ment . 


k 


I 


(t. N o /T) exp(- | y* |/1.253 a ± A ± ) 


(64) 


i = 1 

A^, and a are evaluated separately for each mission seg- 


Dicrete Wave Phenomena . Although planetary scale waves and tides have hor- 
izontal wave lengths much larger than the gravity wave components, their verti- 
cal wavelengths may be comparable. For horizontal flight the presence of these 
other variations would not be important. However, in those cases in which the 
trajectory angle 9 is small but non-zero the effects of the presence of these 
vertical variations can be important. It is suggested that, rather than use a 
discrete wavelength in a spectral approach, it would be more practical to con- 
sider only what mean square variation along the trajectory will occur. Consider 
a simple model in which some phenomenon f (e.g. a wind or thermodynamic variat- 
ion due to tides or planetary scale waves) has a height and horizontal variation 
at a specific time given by 


f (x, z) « A sin(27rx/L^) sin (2 ttz/L^) (65) 

where x is along the trajectory direction, and are the horizontal (along x) 
and vertical wavelengths respectively, A is the amplitude, and the phase was se- 
lected arbitrarily. The mean square differences in the parameter f in going 
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from x to x + r- and z to z + r~ can be determined by 


2 1 f L 3 f L i 

<Af (r x , r 3 )> - — - [f (x + r z + ^ 

13' ' 


- f(x, z) ] dxdz 


( 66 ) 


If equation (65) is substituted into ( 66 ) the result is 

<Af 2 (r 1 , r 3 )> = (A 2 /2) [1 - cos(2ir i^/l^) cos(2tt r 3 /L 3 )] (67) 

In order to estimate mean square variations in winds, pressure, density and 
temperature due to tides and planetary scale waves in the 45 - 65 km region the 
magnitude values listed in Table 10 may be taken as representative. The plane- 
tary wave magnitudes shown in Table 10 are from Figure 2 and the tidal magni- 
tudes are from the analysis results of Section 5. The most appropriate verti- 
cal wavelengths L^, in the 45 - 65 km height range are: approximately infin- 

ite for the semidiurnal tide and the mid-and high latitude diurnal tide, and 
about 30 km for the low latitude diurnal tide. For the planetary scale waves 
in the 45 - 65 km height range - 37 km may be used. Horizontal wavelengths 
for both the tides and planetary waves would be a few thousand kilometers, de- 
pending on the trajectory angle and the latitude. 

The time structure function data shown in Figure 1 indicates fairly large 
magnitude gravity waves with a discrete period of 3 hours. The effects of 
this gravity wave mode may be accounted for by the same method just described 
except that the time variation must also be taken into account, since appre- 
ciable variation of this gravity wave component could occur during the time 
span of re-entry. Equation (67) may be generalized to accout for a time 
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TABLE 10. MAGNITUDES OF TIDAL AND PLANETARY WAVE 
VARIATION IN THE 45 - 65 KM HEIGHT REGION 
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difference x by the relation 


<Af (r 1> r 3 , x)> = 

(A 2 /2) [1 - cos(27rr 1 /L 3 > cos ^nr^/L^) cos(2ttt/P)] 


( 68 ) 


where the period P = 3 hours in this case. The amplitudes of the 3 hour grav- 
ity wave components in the 45 - 65 km height range are estimated from Figure 1 

to be A = 4.2%, A = 4.7%, A_ = 3.2%, A * A = 5.7 m/s. The horizontal scale 
p * p * T u v 

2 2 1/2 

L, is estimated to be L. = (A + A ) P = 90 km and the vertical scale L_ 

1 1 u v 3 

is probably in the range of 5 to 10 km. 

A Re-entry Trajectory Model . Consider a model similar to that of Campbell v "'^ / , 
which divides the trajectory x - z plane into discrete coordinate points x^, z^ 
where x^ = x q + i Ax and z^ = z^ + jAz. We consider that a parameter, say den- 
sity, evaluated at these points would be given by 


(55) 


p . . - p . . + p . . + r . . 

ij K i 3 i 3 


(69) 


where p is the mean density appropriate to the height, latitude, and month or 
season as determined by some reference or standard atmosphere. The component 

I 

p ^ is a discrete wave component and may be evaluated by 


p.. = p + A {[1 - cos(2iriAx/L ) cos(27rjAz/L«) ]/2} 

1J OO P 1 1 -5 


1/2 


(70) 


as from equation (67), or (68) for the 3 hour gravity wave component. repre- 

sents the continuous spectrum gravity wave component and might be determined as 
a set of random numbers selected from a gaussian distribution with the standard 
deviation determined from Table 8. However, the values at neighboring points 
must be correlated. The set of values can be made to have explicit corretion by 
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generating 

relation 


values from surrounding, previously generated values by the 


F i + 1, j + 1 


i, i + 1 


(71) 


1 X X, J 


where e is an uncorrelated gaussian random variable with standard deviation given by 
Table 8 and a, 8, and y are determined from the known correlation between the 
adjacent points (given by equation (37) for scalar quantities or equation (40) , 

(56) 

(41), or (42) for velocity components) by methods developed by Hicks and Justus . 
The spacecraft is then "flown" through an ensemble of sets of such density fields 
and the variation of trajectory parameters can be computed. 
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APPENDIX A. 


ERROR ESTIMATE IN STANDARD DEVIATION 
FROM MEASURED FOURTH MOMENTS 


Consider a set of values X^, i * 1 to n which have mean zero (x = <x> = 0) 


and 


variance 


2 2 2 

a (<x > = a ), We now wish to compute an estimate of 


the error e in the measued value of a. The variance of the set of data values 

o 

2 2 
is v * <x > and the individual contributions to the variance are v. * x. . The 

l i 

standard deviation S of the variance values is given by 


<(v i - v) 2 > 


2 -2 
<V > - V 


(A-l) 


and the error e in the mean variance estimate can be taken as 
v 


(<v 2 > - v 2 )/(N - 1) 


(A-2) 


where N is the number of independent values in the total of N values of x. Since 
2 2 4 4 

v * a and <v > * <x > » 8 o , where 8 is the kurtosis of the x distribution 

4 4 

(i.e. 8 = <x >/o ), then (A-2) becomes 


s v 2 = q 4 (8 - 1)/(N - 1) 


(A-3) 


Now the error e in the mean a value is related to e by 
a v J 


(a + e 0 )‘ 


v + e 
— v 


(A-4) 


but, from expansion of the left side of (A-4) 


(a + e^) 2 - o 2 + 2 o 


(A-5) 


where the second order term in e has been neglected. Since v = a , then it is 

a 
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apparent from combination of (A-4), (A-5) , and (A-3) that 

e o = zjlo = (a/2) [(6 - 1)/(N - 1)] 1/2 (A-6) 

The use of the divisor N - 1 in (A-2) means that the error is of the nat- 
ure of an error of the mean. Thus from (A-6) gives an error which would 

represent the range of deviation of the mean values obtained from comparable 
data sets each made up of N independent values. 
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